Momentum distribution of a Fermi gas of atoms in the BCS-BEC crossover 
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We observe dramatic changes in the atomic momentum distribution of a Fermi gas in the region of 
the BCS-BEC crossover. We study the shape of the momentum distribution and the kinetic energy 
as a function of interaction strength. The momentum distributions are compared to a mean-field 
crossover theory, and the kinetic energy is compared to theories for the two weakly interacting limits. 
The temperature dependence of the distribution is also presented. 



Recent years have seen the emergence of an intriguing 
Fermi system achieved with ultracold gases of 40 K or 6 Li 
atoms. With these systems it is possible to widely tune 
the interatomic interaction strength, represented dimcn- 
sionlcssly as kpa, where kp is the Fermi wavevector and 
a is the scattering length. Of particular interest is the 
strongly interacting regime (—1 < 1/kpa < 1) where a 
crossover between BCS theory of superconductivity and 
Bose-Einstein condensation (BEC) of molecules occurs 
IH Experiments have shown that these Fermi sys- 
tems cross a phase transition as a function of temperature 
and display features of the BCS-BEC crossover such as 
a pairing gap, A, on order of the Fermi energy, Ep. Ex- 
perimental probes have been numerous and include stud- 
ies of molecule formation 0, 00, 01 > thermodynamic 
properties Hoi Jill Jl2L lis! condensate formation 
in the crossover |l5l Il6j . collective excitations fl7l Hif , 
single-particle excitations 0, |2(| , and vortices |2lj . 

One classic phenomenon associated with pairing in a 
Fermi system that has yet to be explored fully in atomic 
systems is a broadening of the Fermi surface in momen- 
tum space (see for example H2)- Figure [I] (inset) shows 
the expected momentum distribution of a homogeneous, 
zero temperature (T = 0) Fermi system. In the BCS 
limit (1/kpa —> —oo) the amount of broadening is small 
and associated with A. As the interaction increases this 
effect grows until at unitarity (1/kpa = 0) the effect is on 
order of Ep, and in the BEC limit (1/kpa — » oo) the mo- 
mentum distribution becomes the square of the Fourier 
transform of the molecule wavefunction. 

In this Letter we present a measurement of the atom 
momentum distribution of a trapped Fermi gas in the 
BCS-BEC crossover regime. To perform these experi- 
ments we prepare a 40 K gas near a scattering resonance 
known as a Feshbach resonance, where we can tune the 
s-wave interaction between fermions by varying the mag- 
netic field B 0, . To probe the system we use the 
standard technique of time-of-flight expansion followed 
by absorption imaging |24). To measure the momentum 
distribution of atoms the gas must expand freely without 
any interatomic interactions; to achieve this we quickly 
change the scattering length to zero for the expansion. 
Bourdel et al. pioneered this type of measurement using 



a gas of 6 Li atoms at T/Tp « 0.6, where Tp is the Fermi 
temperature |l ll| . Here we report measurements down to 
T/Tp w 0.1, where pairing becomes a significant effect 
and condensates have been observed [l5l E| . 

To understand what we expect for our trapped atomic 
system, we can predict the atomic momentum distribu- 
tion using a local density approximation and the results 
for the homogeneous case. In the trapped gas case, in 
addition to the local broadening of the momentum dis- 
tribution due to pairing, attractive interactions compress 
the density profile and thereby enlarge the overall mo- 
mentum distribution. Figure ^ shows an integrated col- 
umn density from the result of a mean-field calculation 
at T = as described in Ref. 0. 

In our experimental setup we create an ultracold K 
gas using previously described cooling techniques 0,H3|- 
The gas is prepared in a nearly equal mixture of the spin- 
states |/, to/) = |9/2, —9/2) and |9/2, —7/2), where / 
is the total spin and m/ the spin-projection quantum 
number. The final ultracold gas is held in an optical 
dipole trap formed at the intersection of two gaussian 
laser beams. One beam is oriented parallel to the force 
of gravity (y) with a waist of w y =200 /im and the second 
beam is perpendicular to the first (z) and has w z =15 /im. 

We first measured the atomic momentum distribution 
with our lowest temperature Fermi gas. We start with a 
weakly interacting gas at T « 0.12 Tp in a trap with a 
radial frequency of v r = 280 Hz and an aspect ratio of 
v z /v r = 0.071. We then adiabatically increase the inter- 
action strength by ramping the magnetic field at a rate 
of (6.5 ms/G) _1 to near a Feshbach resonance located 
at 202.10 ± 0.07 G 0. After a delay of 1 ms, both 
dipole trap beams are switched off and simultaneously a 
magnetic-field ramp to a m (B — 209.6 G) at a rate of 
(2 /is/G) -1 is initiated. The rate of this magnetic-field 
ramp is designed to be fast compared to typical many- 
body timescales as determined by -gj = 90 fis. The cloud 
is allowed to freely expand for 12.2 ms, and then an ab- 
sorption image is taken. The imaging beam propagates 
along z and selectively probes the 1 9/2, —9/2) state fioj . 

Samples of these absorption images, azimuthally aver- 
aged, are shown in Fig. [2 for various values of 1/kpa, 
where the superscript indicates a quantity that was 
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FIG. 1: (color online) Theoretical column integrated momen- 
tum distributions of a trapped Fermi gas n(k) calculated from 
a mean-field theory |2fi||. N is the total particle number, k° F 
is the non-interacting Fermi wavevector. The normalization 
is given by 2ir J n(k)kdk = N, and we note that the nor- 
malization is strongly affected by the tail of the distribution. 
The lines correspond to a = (blue), l/k F a = —0.66 (green), 
1/kpa = (red), and l/k° F a = 0.59 (black), (inset) Corre- 
sponding distributions for a homogeneous system. 



FIG. 2: (color online) Experimental, azimuthally averaged, 
momentum distributions of a trapped Fermi gas normalized 
such that the area under the curves is the same as in Fig. 
The curves correspond to \/k%a = —71 (blue), -0.66 (green), 
(red), and 0.59 (black). Error bars represent the standard 
deviation of the mean of averaged pixels, (inset) Curves for 
l/k° F a = -71 (top) and (bottom) weighted by k 3 . The lines 
are the results of a fit to Eqn. 



measured in the weakly interacting regime. We observe a 
dramatic change in the distribution as predicted in Fig. 
n Some precautions need to be taken in quantitative 
comparison of Figs. ^ and [21 First, the magnetic-field 
ramp to the Feshbach resonance, while adiabatic with 
respect to most time scales, is not fully adiabatic with 
respect to the axial trap period. Second, in the exper- 
iment an adiabatic field ramp keeps the entropy of the 
gas, not T/Tp, constant. However, we expect the result- 
ing change in T/Tp to have a minimal effect on the dis- 
tribution for l/k° F a < @j. Third, the theory assumes 
T = and does not include the Hartree term, thus un- 
derestimating the broadening on the BCS side compared 
to a full theory |29| . 

It is natural now to consider extracting the kinetic en- 
ergy from the momentum distribution. While the mo- 
mentum distribution should be universal for small mo- 
menta, for large momenta it is influenced by details of 
the interatomic scattering potential. In the extreme case 
of a delta potential, which we used for the calculation 
in Fig. 1, the momentum distribution has a tail with 
a I/fc 4 dependence, giving rise to a divergence of the 
kinetic energy. In the experiment we avoid a depen- 
dence of the measured kinetic energy on details of the 
interatomic potential because our magnetic-field ramp is 
never fast enough to access features on order of the in- 
teraction length of the Van der Waals potential, ro ~ 60 
ao for 40 K [30(. Thus, the results presented in this Letter 



represent a universal quantity, independent of the details 
of the interatomic potential. Although universal in this 
sense, the measured kinetic energy is intrinsically depen- 
dent on the dynamics of the magnetic-field ramp, with 
faster ramps corresponding to higher measured energies. 

To obtain the kinetic energy from the experimen- 
tal data exactly we would need to take the second 
moment of the distribution, which is proportional to 
Y.k^OD/Y.kOD, where OD is the optical depth. As 
illustrated in Fig. [21 (inset) this is difficult due to the 
decreased signal-to-noise ratio for large k. Thus, our ap- 
proach will be to apply a 2D surface fit to the image and 
extract an energy from the fitted function. In the limit 
of weak interactions the appropriate function is that for 
an ideal, harmonically trapped Fermi gas. This is 

OD(x,y) =pfc. g2 (-Ce^~^)/. g2 (-C) (1) 

00 k 

where g n (x) = w> ( ' 1S the fugacity, a^. y are pro- 

k=l 

portional to the Fermi gas temperature, and pk is the 
maximum OD. Assuming isotropic expansion in all three 
dimensions the kinetic energy per particle is given by 

_ 3 ma x a y g 4 (-() 
2 t 2 g 3 {-() 

where m is the mass of 40 K and t is the expansion time. 
Empirically, we find Eqn. ^ fits reasonably well to data 
throughout the crossover, as illustrated in Fig. [21 (inset). 
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Figure [3] shows the result of extracting E^ n as a func- 
tion of l/k F a; we see that Ekm more than doubles be- 
tween the non-interacting regime and unitarity. We have 
checked that heating and loss due to inelastic processes 
are negligible up to l/kpa ~ 0. To do this we performed 
an experiment in which we adiabatically approach the 
Feshbach resonance at rate of (6 ms/G) -1 , wait 1 ms, 
and then ramp back at the same slow rate to the weakly 
interacting regime. If we start with a cloud initially at 
T/Tp = 0.10, T/Tp upon return increases by less than 
10% for a ramp to l/kpa = (yet by 80% for a ramp to 
l/k° F a = 0.5). 

Using the fitting function of Eqn. \l\we can also extract 
information about the shape of the distribution through 
the parameter £. Since £ can range from -1 to oo it is 
convenient to plot the quantity ln(l + Q (Fig. 0}. We 
find that the shape evolves smoothly from that of an 
ideal Fermi gas at T/Tp ~ 0.1 in the weakly interacting 
regime, to a gaussian near unitarity, and to a shape more 
peaked than a gaussian in the BEC regime. These quali- 
tative features are predicted by the mean-field calculation 
of the distributions in Fig. 

As mentioned earlier Ej~i n of a trapped gas is affected 
both by the broadening due to pairing (Fig. (inset)) 
and by changes in the trapped gas density profile. In 
the BCS limit, the broadening due to pairing scales with 
e -7r/2fc F |o| anc j j g tjyug exponentially small compared to 
density profile changes, which scale linearly with /cp|a|. 
In this limit we can calculate Ekin lE\j„ using a mean- 
field calculation in the normal state [31| to find, to lowest 
order in k° F \a\, E km /E° km = ^k%\a\ + 1. We plot this 
result in Fig. [3] (inset) and find good agreement for the 
weakly interacting regime (l/k F \a\ > 1). In the crossover 
regime where the pairs are more tightly bound, pairing 
provides a significant contribution to the change in the 
momentum distribution. At unitarity a full Monte Carlo 
calculation predicts the radius of the Fermi gas density 
profile to become (1 + /3) 1 / 4 i? = 0.81i?o, where Rq is 
the Thomas-Fermi radius of a non-interacting Fermi gas 
|2f)j . Just this rescaling would result in Eki n /E® in = 1.54 
(green bar in Fig. Thus, at unitarity, pairing effects 
on the momentum distribution must account for a large 
fraction of the measured value of Eki n /E® in = 2.3 ± 0.3 
(Fig. 0) and all of the observed change in distribution 
shape (Fig. EJ. 

In the BEC limit we expect the measured energy to 
be that of an isolated diatomic molecule after dissoci- 
ation by the magnetic-field ramp. Provided the scat- 
tering length associated with the initial molecular state, 
a(t = 0), is much larger than rg ss 60 ao, the wave func- 
tion for the molecule is given by tp = e -r / a '* =0 '/r where r 
is the intcrnuclcar separation. We can calculate the mea- 
sured energy from the solution of the Schrodinger equa- 
tion with a time-dependent boundary condition on the 
two-particle wavefunction dl °s( r V') _ — 1 wnere 
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FIG. 3: (color online) The points (•) show the measured en- 
ergy E kin normalized to E kin — 0.25 kb /xK (E kin = |-Ef 
for a harmonically trapped gas at T = 0). The red line is 
the expected energy ratio from a calculation only valid in the 
weakly interacting regime (l/kpa < —1). In the strongly in- 
teracting regime pairing due to many-body effects strongly 
increases Ekin- The green bar represents the expected value 
of E^iri. I E® in at unitarity just due to density profile rescaling 
(see text). In the molecule limit (l/kpa > 1) we calculate the 
expected energy for an isolated molecule (blue line), (inset) 
A focus on the weakly interacting regime with the same axes 
definitions as the main graph. The dashed lines show the un- 
certainty in the calculation based upon the uncertainty in the 
Feshbach resonance parameters |15ll23|l . 



a(t) is the scattering length fixed by the magnetic field 
at time t. In Fig. [3|we show the result of this calculation 
for a pure gas of molecules and a (2 /xs/G) -1 ramp rate. 
We find reasonable agreement considering that there is a 
large systematic uncertainty in the theory curve due to 
the experimental uncertainty in the magnetic-field ramp 
rate and that this two-body theory should match the data 
only in the BEC limit (l/k F a > 1). 

We have also studied the dependence of the momentum 
distribution on (T/Tp)° . To vary the temperature of our 
gas, we recompress the optical dipole trap after evapora- 
tion and parametrically heat the cloud [8j. The experi- 
mental sequence for measuring the momentum distribu- 
tion is the same as above except the ramp rate to a = 
for expansion was ~ (8 /xs/G) -1 . Figure [S] shows the 
measured kinetic energy change AEkin = Ekin — E° in . 
On the BEC side of the resonance (l/kpa > 0), AE^m 
decreases dramatically with (T/Tp)°. Because AE^m 
should be proportional to the molecule fraction, this re- 
sult is closely related to the recent observation that the 
molecule conversion efficiency scales with T/Tp 0. In 
the strongly interacting regime we also observe a decrease 
in AE^in with increasing (T/Tp)° . Here the temperature 
dependence of AEkm is consistent with the expectation 
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FIG. 4: Shape of the momentum distribution as described 
through the parameter £ (Eqn. 0. ln(l + £) = corresponds 
to a gaussian distribution, and for an ideal Fermi gas ln(l + 
Q- 1 ?s T/Tf in the limit of low T/T F . The dashed lines show 
the uncertainty in the Feshbach resonance position [TrI |. 
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FIG. 5: (color online) Temperature dependence of AEkin = 
Ekin — -Efcin normalized to the Fermi energy at a = 0. (T/Tf)° 
is the temperature of the non-interacting gas |3^ |. For the 
coldest dataset (black) the peak density, for atoms in one of 
the two spin states, in the weakly interacting regime is Up k = 
1.4 x 10 13 cm~ 3 and E% = 0.56 /uK. For the hottest dataset 
(magenta) n^ k decreases to 6 x 10 12 cm" " 3 and E% = 0.79 yK. 



that the changes in the kinetic energy are caused by pair- 
ing and not coherence ^M^SE3- 

We have found that the momentum distribution of a 
Fermi gas provides a wealth of information on crossover 
physics. The measurement is a probe of pairing in the 
strongly interacting regime, and it provides a univer- 
sal thermodynamic quantity that is complementary to 
reviously measured energies in the BCS-BEC crossover 
Hol fill H^ . Il3| . This work also provides a starting 
point for future experiments seeking to probe pair cor- 
relations in the momentum distribution using measure- 
ments of atom shot noise (34). 
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